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i-j— ♦ ' Abstract - Experiments and numerical simulations reveal that in the forward cascade regime, 

I ' the energy spectrum of two-dimensional turbulence with Ekman friction deviates from Kraichnan's 

^ . prediction of fc^^ power spectrum. In this letter we explain this observation using an analytic model 

t^_H ' based on variable enstrophy flux arising due to Ekman friction. We derive an expression for the 

^j^ , enstrophy flux which exhibits a logarithmic dependence in the inertial range for the Ekman-friction 

O ' dominated flows. The energy spectrum obtained using this enstrophy flux shows a power law scaling 

' ly-j ' for large Reynolds number and small Ekman friction, but has an exponential behaviour for large 

^~^, Ekman friction and relatively small Reynolds number. 



s^j^ \ Physics of turbulent flow is quite complex. One of the important and generic features of 

\^ i turbulent flow in three dimensions (3D) is a constant energy flux from large length scales 

^S| \ to small length scales. The wavenumbers exhibiting constant energy flux have k~^'^ energy 

VO ' spectrum [l]. Two-dimensional (2D) fluid turbulence however has significantly different be- 

fT^ I havior. Kraichnan [2] showed that in 2D turbulence, the low wavenumber modes (below the 

^^ ■ forcing wavenumber) exhibit inverse energy cascade, while the large wavenumber modes have 

Cn I forward enstrophy (square of the vertical vorticity) cascade. These two regimes have k~^'^ 

and k~^ energy spectra respectively. The above features have been observed in numerical 
simulations 0^6], and in experiments involving electromagnetically driven flows [Tll^ and 
soap films [QjIlO]. The atmospheric data indicates that some features of the atmospheric 
J, ^ turbulence is two-dimensional [TT|[12]. 

^_' The fluid flow in a two dimensional surface is also affected by the drag from its envi- 

ronment as demonstrated by several experiments |13H17j and numerical simulations J18ill9j . 
This kind of friction is also referred to as "Ekman friction". Belmonte et al. [13] per- 
formed experiments on freely suspended soap films, and reported that the energy spectrum 
is steeper than k~^ . Boffetta et al. [TTj observed similar steepening of the spectrum in their 
experiment on a thin layer of ionic fluids, which is driven electromagnetically and is moving 
under a layer of fresh water. Here, the drag on the ionic fluid is induced by the shear layer 
induced by fresh water [l4l[T7] . Nam et al. [18] and Boffetta et al. [1^ performed numerical 
simulations with Ekman friction and reported the spectral exponent to be larger than 3. 

The aforementioned steepening of the energy spectrum in two-dimensional flow has been 
attributed to Ekman friction, which is modeled as — au, where a is a constant, and u is 
the velocity field. Nam et al. [TH] and Boffetta et al. [T7] attempted to derive the new 
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spectral indices using an analogy with the dynamics of a scalar in a turbulent fluid. They 
postulated that a fluid blob has a finite lifetime r = a'~^, and it is passively advected. The 
turbulent motion leads to a stretch of a fluid blob with a mean rate given by the Lyapunov 
exponent A. Using the above assumption, the incompressibility condition, and the mean field 
approximation, Boffctta et al. [T7] predicted that the energy spectrum E{k) ^ fc""^"^"/^. 
Boffetta et al. ^9] also studied the intermittency effects of Ekman friction, and showed 
that the small-scale statistics of the vorticity fluctuations is related to the passive scalar 
transported by the velocity field. Perleker and Pandit |20j performed detailed numerical 
simulation and studied the effects of Ekman friction on the structure function. They showed 
that the velocity structure functions display simple scaling. The reader is also referred to a 
review article by Kellay and Goldburg [16]. 

The arguments of Boffetta et al. [17] involves several assumptions and relatively com- 
plex computation of the Lyapunov exponent. In the present letter we compute the energy 
spectrum using a model based on a variable enstrophy flux. Unlike the viscous force, which 
affects the dissipation range, Ekman friction is active at all scales. As a result, the enstro- 
phy flux decreases significantly in the inertial range itself. This decrease in the flux leads 
to a steepening of the energy and enstrophy spectra. In the present letter, we will derive 
an expression for the variable enstrophy flux in terms of the dissipative parameters a, the 
kinematic viscosity v, and other turbulence parameters. We will show that the variable 
enstrophy flux yields energy spectrum steeper than k~'^. 

The dynamical equations for the vertical vorticity in a two-dimensional fluid flow with 
Ekman friction are 

-^ + u-Vu} ^ -au} + iy\/'^uj + fit), (1) 

ot 

and the incompressibility constraint V • u = 0, which implies that the density of the fluid 
is a constant. In the above equation, u is the velocity field, a; = (V x u)^ is the vertical 
component of the vorticity field, and v is the kinematic viscosity of the fiuid. An external 
force f{t) is applied to maintain a steady-state. 

The enstrophy fiux n(fc) is defined as the total enstrophy transferred from the modes 
inside the wavenumber sphere of radius k to the modes outside the sphere. Since the above 
flux is defined in the wavenumber space, we rewrite Eq. ([T]) in the Fourier space as 

^^ = -ik, Y, zi,(q)c^(k - q) - acoik) - i^k^^) + /(k), (2) 

q 

where k is the wavenumber. The corresponding enstrophy evolution equation is 

dZ(k) 



dt 



T{k)-2aZ{k)-2vk^Z{k)+F{k), (3) 



where Z{k) = |a;(k)p/2 is the one-dimensional enstrophy spectrum, T(fc) is the enstrophy 
transfer term arising due to nonlinearity, and F(k) is the enstrophy supply rate due to the 
forcing [T]. It has been analytically shown that the total energy E = 'Y^ |u(k)p/2 and the 
total enstrophy Z ^ ^ |a;(k)p/2 are conserved in the inviscid limit (z^ = and a = 0) 
without any external forcing [l]|2] . The dual energy spectrum described in the introduction 
is a consequence of simultaneous conservation of energy and enstrophy [2] . Two-dimensional 
turbulent flow with a = and v ^i- Q exhibits k~^^^ energy spectrum for k < kf, and k^^ 
energy spectrum for k > kf, where kf is the center of the forcing- wavenumber band. The 
scenario changes significantly when a 7^ 0, as observed in earlier experiments and numerical 
simulations. 

The forcing wavenumber in a typical experiment on Ekman friction is of the order of 
the box size. Consequently, experiments report only forward enstrophy cascade and energy 
spectrum as fc^° with a > 3. Since \oj{k)\ = fc|u(k)|, the enstrophy spectrum would vary as 
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^-a+2_ Pope [1] prescribed a function for the energy spectrum for three-dimensional fluid 
turbulence. We adapt Pope's function for the enstrophy spectrum Z{k) as 

Zik) - C(n(fc))2/3fc-VL(fci)A(fcry), (4) 

where C ~ 1.4 is a constant (equivalent to Kolmogorov's constant for 3D fluid turbu- 
lence) [2T1[22], n(fc) is the enstrophy flux emanating from the wavenumber sphere of radius 
fc, and /i(fcL),/^(fc?7) specify the components of the forcing-scale and dissipative-scale en- 
strophy spectra respectively. These functions have been described by Pope [T] as 

f„(kv) = e,Lp[-ll{l(k,i}' + 4f"-c„}], (6) 

where L is the box size, and cl, c^,Poj /? SlTC constants, which are determined by matching the 
function of Eq. Q with experimental observations. In the absence of any clear prescription 
for these constants, we take Cl ~ 6.78, c,, w 0.40, /3 ~ 5.2 and po = 2, used by Pope [I] for 
3D turbulence. 

In Kraichnan's phenomenology for the two-dimensional fluid turbulence without Ekman 
friction, the enstrophy flux n(fc) is a constant in the inertial range. This is due to the fact 
that in this regime, the local dissipation rate {2vk^Z{k)) is negligible, and the forcing is 
absent. However, in the presence of Ekman friction, n(fc) decreases in the inertial range 
itself due to the 2aZ{k) term. The enstrophy fluxes n(fc -I- dk) differs from n(fc) by the rate 
of enstrophy loss in the wavenumber shell (fc, k + dk), i.e., 

n(fc + dk) - U{k) = - {2iyk^ + 2a} Z{k)dk, (7) 

or 

S) ^ _ {2^fc2 + 2a] C{n{k))^'^k-'f,{k^). (8) 

In the present letter, we focus on the inertial and dissipative ranges. Therefore, we have 
ignored the variation of the forcing-scale spectrum /^(fcL), i.e., /^(fci) « 1. 
Equation ([8]) can be easily integrated, which yields 

Ii{k)V^ . 2C V _, . 2aC 



Hn 



Q 1/3 o^i^^'^) r73^2(fc??) 

= l-^Ci/i(fcr;)-^/2(fcr;), (9) 

where Ho is the maximum value of the enstrophy flux at A: = k\ (which is also the lower 
limit for the integration), C\ = i^/(nQ r^) is a constant, and 



is the Kolmogorov length for two-dimensional turbulence. For all our calculations, we also 
take the constant C to be 1.4 [21 ] l22j . The integrals /i and I2 are 



hik-n) = / dk'k'f,,ik') (11) 

'J kiT] 

hikil) = / dk'k'-^fr,{k') (12) 

J k\r} 
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where the lower hmit for the wavenumber ki is chosen to be 6 x 27r/L [T] . Pope [T] suggests 
that the incrtial range starts around this wavenumber. Using Eq. (|10p and Hg = U^ /L^, 
we can deduce that kirj = 6{2Tr/ L)ri = l27r/^/CiRe. In our model the enstrophy flux peaks 
at k = ki and decreases for higher values of k. Equation ([9]) shows that the enstrophy flux 
depends critically on the nondimensional parameter a' = 2aC / {3IIq ). 

The form of Eq. ([9]) appears quite complex, but it can easily approximated in the inertial 
range in the t^ — > limit. Under this limit, /i(fc) w and hik) ~ a'log(fc/fci). Therefore, 
the approximate enstrophy flux would be 

ffl«[l-a'log(fc/fci)f, (13) 

tio 

which shows that for nonzero a', the enstrophy flux decreases logarithmically with the in- 
crease of wavenumber. We will show later that for turbulent flows, the above approximate 
function matches quite well with that computed using Eq. ^ in the inertial range. Ad- 
ditionally, the wavenumbers satisfying the condition a'log(fc/fci) <C 1 follow a power law 
behavior: 



Un 



g-da log(fc/tei) _ (fc/fci)-^" (14) 



Once we have computed the fc-depcndent enstrophy flux, the enstrophy spectrum can be 
immediately derived as 



cn^/'fc-i/„(M 



n(fc) 



2/3 



if fc > fci 



^(fc)=<( 1' -"n-'/^LnoJ ' "— - (15) 

CIIq' k ^/^{kL), otherwise. 



while the energy spectrum varies as 



n(fc) 



if fc > fci 



2/3 

E{k) = <; ^'-^S J''^'" L^J ' " (16) 





Cnl^^k-^fL{kL), otherwise. 



For a narrow range of the inertial range with a' \og{k/ki) <C 1, the aforementioned equations 

yield 

£;(/fc) - Cno/\"-3-2a; (^7) 

We can also compute the enstrophy dissipation rates due to the viscous force and Ekman 
friction using Dy{k) = J^ 2vk^Z{k)dk and Da{k) = J^ 2a Z{k)dk respectively. Another 
important quantity is the ratio of the nonlinear term and Ekman friction, which can be 
estimated using Rbq = U/{aL), termed as "Reynolds number based on Ekman friction". 
We will show later that Rcq plays an important role in determining the structure of the 
energy spectrum. 

We use the aforementioned model to study the nature of energy spectrum as a function of 
the Ekman friction parameter a. We take two cases: large Reynolds number, and relatively 
small Reynolds number. For the former case, we take Re = 10^ and a' = 0, 0.3, and 0.5, 
and compute the normalized enstrophy flux n(fc77)/no and energy spectrum E (ki]) / E (kirj) 
using equations ([9]) and (J16p respectively. These functions are exhibited in Fig. [1] using 
red {a' ~ 0), green {a' — 0.3), and blue {a' ~ 0.5) curves. For a' = 0, the incrtial range 
enstrophy flux is nearly constant and E{k) ~ k~^, consistent with Kraichnan's predictions, 
more so a validation of our model. Note however that n(fc) must vanish for wavenumbers 
k ^ rj^^. We use this condition to determine the constant Ci of Eq. ©, which yields 
Ci = 7.17. 

With the increase of a' or Ekman friction, the entropy flux becomes fc-dependcnt and 
the energy spectrum is steeper than k^^. For a' = 0.3 and 0.5, the spectral indices are 
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Fig. 1: (a) Plots of energy spectra for Re = 10^ and a' = (red thin curve), 0.3 (green thicker 
curve), and 0.5 (blue thickest curve). The best fit curves fc^" with a = 3,3.8 and 4.2 are plotted 
as dashed cyan lines. The approximate energy spectrum using the approximate enstrophy flux of 
Eq. (|13p is exhibited as chained lines, (b) Plots of the corresponding enstrophy flux (solid lines), the 
approximate enstrophy flux of Eq. (|13p (chained lines), and its power law approximation described 
by Eq. (|14[l (cyan dashed lines). 



3.8 and 4.2 respectively. The best fit curves are shown in Fig. [TJa) using the dashed cyan 
lines. We also plot the approximate enstrophy flux of Eq. ()13|) and the corresponding energy 
spectrum in Fig. [l]as chained lines. In the inertial range, the approximate energy spectrum 
and the approximate enstrophy flux fit rather well with the model predictions. For a very 
narrow range of fcr;, Il{k) ~ k~^°' as shown by the cyan dashed line in Fig.[ljb). This result 
indicates that the the logarithmic dependence of enstrophy flux is valid for a larger range of 
wavenumbers than the power law of Eq. (|14p . This observation is also consistent with the 
fact that the spectral index of the energy spectrum is not exactly equal to 3 + 2a' [Eq. pT]) ] , 
but these two exponents are reasonably close. 

We also compute Reo, for a' = 0.3 and 0.5 and find them to be approximately 3 and 2 
respectively (see Table 1). This result indicates that the nonlinear term u • Vu dominates 
Ekman friction, consistent with the power law behaviour of the energy spectrum. The ratio 
Da/Du is 22 and 100 for a' = 0.3 and 0.5 respectively, indicating that the Ekman friction is 
stronger than the viscous term. The strong wavenumber dependence of the enstrophy flux 
is due to the important role played by Ekman friction. 

We perform our model calculations for the experimental parameters of Boffetta et al. |17j . 
The Reynolds number for these runs ranges from 3950 to 6600, which is relatively small 
compared to the aforementioned test cases with Re = 10^. In Boffetta et al.'s experiment, a 
thin layer of electrolyte solution of water and NaCl contained in a square tank (L = 50 cm) 
is stirred electromagnetically using four triangular shaped electromagnets with alternating 
signs. A layer of fresh water was kept above the electrolyte solution. The thickness of the 
ionic fluid was maintained at 0.3 cm, while that of the upper fluid was varied from 0.3 cm 
to 0.7 cm. The ionic fluid experiences a drag force —au. The kinematic viscosity of the 
electrolyte is same as that of water {v sa 0.01 cm^/sec). For the three sets of experiments 
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performed, Boffetta et al. reported the Ekman friction coefficient a = 0.037, 0.059, 0.069, 
Re = u^^sL/v = 6600, 6650, 3950, and tlie spectral indices a = 3.5, 3.8, and 4.0 respectively. 
See Table 1 for more details. 




Fig. 2: Plots of the energy spectra E" (kri) / E" (kirj) computed using our model for the experimental 
parameters of Boffetta et al. [T7] on a log-log scale. The red, green, and blue curves represent the 
Q = 0.037,0.059,0.069 experimental runs respectively (see Table f). The power law function A:^" 
(dashed curves) does not fit with the energy spectrum. 
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Fig. 3: (a) Plots of energy spectra for the experimental parameters of Boffetta et al. [17] with 
a' — 0.037 (red thin curve), 0.059 (green thicker curve), and 0.069 (blue thickest curve) on a semi- 
log scale. See Table 1 for more details. Exponential functions exp{—b{kr] — kir])) (dashed cyan lines) 
with b — 16.7, 18.1, and 16.8 respectively are the best fit curves for the spectra in the non-viscous 
range, (b) The corresponding enstrophy flux for the aforementioned experimental cases. 
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For the three aforementioned experimental sets of Boffetta et al. [17] , the nondimensional 
parameters are a' = 0.046,0.086,0.107 and Re„ = U/{aL) = 0.71,0.45,0.22 respectively, 
indicating that the nonlinear term is somewhat weaker than the Ekman friction. We also 
estimate the maximum enstrophy flux Hg as lu'^j-^^/T w ^rms- Using the aforementioned 
parameters we compute the energy spectra and the normalized enstrophy fluxes for these 
experimental sets using our model, and compare the results. In Fig. [2] we plot the energy 
spectra and their best power-law fits. The figure illustrates that a power law is not a good 
fit for the energy spectra predicted by the model. Therefore, we attempt to fit the energy 
spectra with an exponential function exp(— 6(fc77 — kir])) for k > ki. In Figure |3] we plot 
the normalized energy spectra E(kr])/E(kiri) and normalized enstrophy flux on a semi-log 
scale. The exponential functions (the dashed lines) fit well with the energy spectra in the 
non-viscous range, and the fitting parameters b for the three experimental cases are 16.7, 
18.1, and 16.8 respectively. Thus an exponential function may describe the energy spectrum 
of the Boffetta et al. better than a power-law function, consistent with the observation 
that Ekman friction dominates the nonlinear term in these cases (Rcq. < 1). Also, the 
approximate enstrophy flux computed using Eq. (|13l) does not match with the model flux 
[Eq. ([9])] due to absence of an inertial range, which is due the dominance of Ekman friction 
over the nonlinear term. A cautionary remark however is in order: our estimate of Ho, Re 
using the experimental parameters of Boffetta et al.'s [T7] have significant errors. Hence 
the model predictions are somewhat uncertain, and they need to be verified carefully using 
more refined experiments and numerical simulations. 

In summary, we showed that Ekman friction strongly affects the enstrophy flux in two 
dimensional turbulence. We derived an expression for the variable enstrophy flux in terms 
of Ekman friction parameter a. As a result, the kinetic energy spectrum is steeper than 
k~^. For large Reynolds number flows with Rea = U/{aL) > 1, we expect a power law 
behaviour for the energy spectrum. This is due to the dominance of the nonlinear term over 
Ekman friction. However, Ekman friction is stronger than the nonlinear term for flows with 
Rcq, < 1, hence an exponential energy spectrum may be expected for such flows. 

The method presented here is quite general and could be applied to other situations as 
well. The liquid metal flows in the presence of a strong magnetic field under the quasistatic 
approximation |23) has a dissipation term similar to the Ekman friction. Adaption of the 
present procedure to the liquid metal flow would provide valuable insights into its energy 
spectrum and dynamics. Another important direction for future work could be verification 
of our model using direct numerical simulation. 

I thank Mani Chandra for Python related help and useful comments. I also thank Sagar 
Chakraborty, Prasad Perlekar, Stefano Musacchio, and Anupam Gupta for the discussions. 
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List of Tables 



Table 1: Table depicting our model calculation for the Reynolds number Re = 10^ and a' = 
0,0.3,0.5 (first three rows), and for the experimental parameters of Boffetta et al.'s [17] (next 
three rows). Here a is the coefflcient of the Ekman friction, Umis is the rms velocity, ojims is 
the rms vorticity, Ho is the maximum value of enstrophy flux, a' — aC/(3IlQ ) with C — 1.4, 
Rea = U/{aL) is the ratio of the nonlinear term and the Ekman friction term, a,aB are the 
spectral exponents computed using our model and by Boffetta et al. respectively, b is the coefficient 
of E{k) — exp{—b{kri — kirj)), and Da/Dv is the ratio of the dissipation rates by Ekman friction 
and the viscous force. For Boffetta et al.'s experimental parameters [17], we use Rea ~ Umis/iaL) 
with L = 50 cm, and Ho = cj'^ 



j3 



a 

s-i 


'^rins 

cm/s 


s-1 


Ho 

S-3 


a' 


Re 


Re^ 


a 


as 


b 


Da/D, 


- 


- 


- 


- 





10" 


oo 


3 


- 


- 





- 


- 


- 


- 


0.3 


10*^ 


«3 


3.8 


- 


- 


22 


- 


- 


- 


- 


0.5 


10^ 


w2 


4.2 


- 


- 


100 


0.037 


1.32 


0.75 


0.42 


0.046 


6600 


0.71 


- 


3.5 


16.7 


0.091 


0.059 


1.33 


0.64 


0.26 


0.086 


6650 


0.45 


- 


3.8 


18.1 


0.18 


0.069 


0.79 


0.60 


0.21 


0.107 


3950 


0.22 


- 


4.0 


16.8 


0.16 



p-9 



